Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia.

BACKGROUND: The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia. METHODS: Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data. RESULTS: Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus' populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107-0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Comparative cytogenetics of the ground frogs Eupsophus emiliopugini Formas, 1989 and E. vertebralis Grandison, 1961 (Alsodidae) with comments on their inter- and intraspecific chromosome differentiation.

South American frogs of the genus Eupsophus Fitzinger, 1843 comprise 10 species. Two of them, Eupsophus vertebralis Grandison, 1961 and E. emiliopugini Formas, 1989 belong to the Eupsophus vertebralis group, exhibiting 2n = 28. Fundamental number differences between these species have been described using conventional chromosome staining of few specimens from only two localities. Here, classical techniques (Giemsa, C-banding, CMA3/DAPI banding, and Ag-NOR staining), and fluorescence in situ hybridization (FISH, with telomeric and 28S ribosomal probes), were applied on individuals of both species collected from 15 localities. We corroborate differences in fundamental numbers (FN) between E. vertebralis and E. emiliopugini through Giemsa staining and C-banding (FN = 54 and 56, respectively). No interstitial fluorescent signals, but clearly stained telomeric regions were detected by FISH using telomeric probe over spreads from both species. FISH with 28S rDNA probes and Ag-NOR staining confirmed the active nucleolus organizer regions signal on pair 5 for both species. Nevertheless, one E. emiliopugini individual from the Puyehue locality exhibited 28S ribosomal signals on pairs 4 and 5. Interestingly, only one chromosome of each pair showed Ag-NOR positive signals, showing a nucleolar dominance pattern. Chromosomal rearrangements, rRNA gene dosage control, mobile NORs elements, and/or species hybridization process could be involved in this interpopulation chromosomal variation.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae).

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Mitochondrial genomes of the South American frogs Eupsophus vertebralis and E. emiliopugini (Neobatrachia: Alsodidae) and their phylogenetic relationships.

We report the sequencing and compare the mitochondrial genomes of the South American ground frogs Eupsophus vertebralis and E. emiliopugini and reconstruct phylogenetic relationships among Eupsophus species. These genomes consist of 16,156 and 16,711 bp in length, respectively and contain 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes (tRNA), and partial non-coding D-loop region. Both genomes share 94.5% identity with 879 variable sites. A phylogenetic analysis with other available mitogenomes recovered both species as the sister clade of Alsodes gargola. Sequences from D-loop, CO1, and Cyt b, amplified and sequenced with primers developed from the mitochondrial genomes, allowed us to reconstruct phylogenetic relationships among Eupsophus species. Since our report represents the first mitogenomes for the genus Eupsophus, we expect these data will be valuable for further studies on conservation genetics and on the evolution of Patagonian amphibians.